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Magnetisation dynamics in exchange spring magnets have been studied using simulations of the 
FePt/Fe bilayer system. The FePt hard layer exhibits strong perpendicular magnetocrystalline 
anisotropy, while the soft (Fe) layer has negligible magnetocrystalline anisotropy. The variation of 
the local spin orientation in the Fe layer is determined by the competition of the exchange coupling 
interaction with the hard layer and the magnetostatic energy which favours in-plane magnetisation. 
Dynamics were studied by monitoring the response of the Fe layer magnetisation after the abrupt 
application of a magnetic field which causes the systems to realign via precessional motion. This 
precessional motion allows us to obtain the frequency spectrum and hence examine the dynamical 
magnetisation motion. Since the rotation of the spins in the soft layer does not have a well defined 
magnetic anisotropy, the system does not present the usual frequency field characteristics for a thin 
film. Additionally we obtain multi-peaked resonance spectra for the application of magnetic fields 
perpendicular to the film plane, though we discount the existence of spin wave modes and propose 
that this arises due to variations in the local effective field across the Fe layer. The dynamic response 
is only considered in the Fe layer, with the FePt layer held fixed in the perpendicular orientation. 



I. INTRODUCTION 

Increasing attention has been paid to the so-called 
exchange- spring systems, which couple the hard and soft 
magnetic properties of the two magnetic components via 
the exchange interaction between them. Recent stud- 
ies have indicated that the magnetic energy product can 
be theoretically as high as 120 MGOe pQ. The exchange 
spring system is in general characterised by the reversible 
dynamics in the soft component arising from the compe- 
tition between the magnetic anisotropics in the two mag- 
netic materials. Goto et al. [2] were the first to consider 
such systems based on the assumption that the exchange 
and anisotropy energies are large enough to overcome 
the magnetostatic energy of the system. The specific 
spin configuration for such a system will depend on the 
relative strengths and directions of the anisotropics, as 
well as the exchange interaction between the magnetic 
layers (phases) and the orientation of any applied mag- 
netic field. The spin configuration is then determined 
by the energy minimisation, taking into account the var- 
ious magnetic contributions in the system as a whole. 
Kneller and Hawig [3 applied these principles to the con- 
struction of a new class of permanent magnetic materials 
which are referred to as the exchange springs. These au- 
thors analysed the reversible component of the hysteresis 
loop. Fullerton et al [4] considered layered systems of 
Sm-Co/Fe and Sm-Co/Co, performing simulations based 
on a ID model. Bowden et al. [5] have made exten- 
sions to Goto's ID model for discrete systems and have 
obtained a relation between the bending and exchange 
fields. Bowden et al. PQ have also performed numerical 
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simulations of 1, 2 and 3D systems, using point defects in 
the chain to determine the behaviour of the system and 
predict a collapse of the exchange spring system. The 
model of Asti et al. [6] considers a FePt (hard) mag- 
netic layer with perpendicular anisotropy coupled to a 
Fe (soft) magnetic layer. Among the various consider- 
ations, the authors have studied the equilibrium phase 
diagram as a function of the thicknesses of the two mag- 
netic layers, predicting a transition from a rigid mag- 
net (RM) to an exchange spring (ES). The rigid magnet 
regime concerns the state where the soft magnetic layer 
magnetisation strictly follows that of the hard layer. The 
exchange spring effect is manifest in the local variation 
of the effective field. In the bilayered system that we 
are considering, the FePt layer has dominant perpendicu- 
lar magnetocrystalline anisotropy, while the Fe layer has 
dominant in-plane anisotropy due to shape effects. As 
the Fe layer thickness increases, the local field reflects 
the varying influences of the two contributions and the 
local spin orientation is seen to rotate from the perpen- 
dicular direction at the FePt/Fe. 

Ferromagnetic resonance can provide explicit infor- 
mation regarding interlayer interactions in multilayered 
structures and has been successfully applied to many sys- 
tems [7] and [8]. There are very few papers that deal 
explicitly with FMR in exchange spring systems, so it 
is therefore a new challenge that needs to be studied. 
Grimsditch et al [9 studied experimentally the magnon 
frequency dependence of a Sm-Co/Fe system suing Bril- 
louin scattering and found a good agreement with the 
theory developed for the one-dimensional model. Also, 
Pechan et al [10] investigated the frequency dependence 
of the ferromagnetic modes as a function of the in-plane 
angle. 

Crew and Stamps [11 have simulated a bilayered sys- 
tem of CoPt/Co, to study the resonance frequency spec- 
trum and anisotropy for a particular applied field direc- 
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tion. However, this paper offers a more detailed descrip- 
tion of the dynamical behaviour of these kind of systems, 
contributing to a better understanding. 

Experimentally, FePt/Fe system has been studied by 
Casoli et al. [12], who have performed measurements of 
the hysteretic behaviour, confirming the transition from 
rigid magnet to exchange spring behaviour. Schmool et 
al. [13 have performed ferromagnetic resonance mea- 
surements on the same FePt/Fe interface system and 
have observed significant changes in the FMR spectra 
for RM and ES samples. 

Exchange spring systems offer various potential tech- 
nological applications, for example, in the form of per- 
manent magnets [3], mag-MEMS [6 and MRAM devices 
[T4] . In the current paper we consider the exchange 
spring system which corresponds to the FePt/Fe bilayer 
samples discussed above, in which the magnetocrystalline 
anisotropy of the FePt layer aligns its magnetisation in 
the direction perpendicular to the film plane. 

The paper is organised as follows: in the following sec- 
tion we discuss the basic theory and micromagnetic simu- 



lations used to study the system, in|III]we present results 
and discussions of these simulations and in the final sec- 
tion we will give the conclusions of our study. 



II. THEORETICAL MODEL 



A. Basic theoretical considerations 



and is related to the free energy of the system as de- 
scribed by the relation: 
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The boundary conditions of a magnetic film can have 
a significant influence on the dynamic response of the 
magnetisation to an excitation field or torque [15]. This 
is because these boundary conditions are related to the 
surface and interface anisotropics and can effectively pin 
the magnetic spins at such magnetic discontinuities. For 
example, in a thin magnetic film with a magnetic field 
applied perpendicular to the plane standing spin waves 
can be excited, giving a resonance condition which is de- 
scribed by: 



ujn = 7[H -H k -M + Dkl] 



(5) 



where k n is the spin wave wavevector of mode n. Restric- 
tions imposed by the boundary conditions on the allowed 
values of the wavevectors can be expressed as [T6] : 



k (pi +P2) 



Pi P2 
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Pi and p2 are the pinning parameters which depend on 
the surface or interface energies and hence will be related 
to surface and interface anisotropics. 



The free energy of a magnetic system can be expressed 
in the most general form as the sum of various contri- 
butions: exchange energy, E exi demagnetising energy, 
Edemj magnetocrystalline anisotropy energy, Ek and the 
Zeeman energy, Eq. This can then be expressed as: 

E tota i = E ex + E dem + E + E K • (1) 

The equilibrium condition of the magnetic system can 
frequently then be evaluated by minimising this energy 
with respect to the orientation of the magnetisation for 
example. The dynamics of coherent magnetisation be- 
haviour is typically described using the Landau - Lifshitz 
- Gilbert (LLG) equation and takes the form: 




This equation describes the temporal evolution of the 
direction of the magnetisation vector due to the preces- 
sional motion induced by the magnetic torque associated 
with the effective applied field, H e ff and its subsequent 
relaxation (viscous damping) which is controlled by the 
magnitude of the Gilbert damping parameter, a. The 
effective magnetic field will have various contributions: 

~^eff ~ + ^0 + ^ex ~ ^ dem , (3) 



B. Simulations 

In this work we concentrated on the study of the static 
and dynamic properties of a bilayer system where, for 
the most part, the interface has a rigid interface condi- 
tion. This means that the interface spins at the Fe/FePt 
boundary are fixed in the perpendicular direction (which 
we take as the y-axis of the coordinate system). The 
dynamic response comes entirely from the Fe layer for 
our simulations (which we introduce in light of experi- 
mental results, which show that this is active in the fre- 
quency range studied [13 ]). The simulations were per- 
formed by constructing a bilayer of FePt, with perpen- 
dicular anisotropy (in the Llo phase [HJ[T7]), coupled to 
an Fe layer using the OOMMF (Object Oriented Micro- 
magnetic Framework) [18] software. The Fe film thick- 
ness was variable, in general, but for the dynamics we 
have maintained this constant, once we have predefined 
the thickness necessary to obtain an exchange spring sys- 
tems, which corresponds to 25 x 10 -10 m. We note that 
the lateral dimensions were arbitrarily chosen and are 
much larger than the film thickness. The Fe layer is con- 
sidered to have shape anisotropy and negligible magne- 
tocrystalline anisotropy (as compared to the FePt layer). 
A Zeeman contribution will arise from the application of 
a magnetic field. 

The precision of the numerical results are governed by 
the micromagnetic cell parameters; (10000, 1, 10000) x 
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10 -10 m. This essentially means that we are taking the 
cell to be one spin in the direction perpendicular to the 
film plane, where we measure the direction with respect 
to the xz plane orientation, as illustrated in fig [I] The 
exchange constant between Fe/Fe and FePt/Fe spins is 
1 x 10 -12 J/m, the uniaxial anisotropy in Fe/Pt layer 
is 4 x 10 6 J/m 3 (in y direction) [19], the damping coef- 
ficient is 0.01 and the time step is 1 x 10 _14 s. Finally 
the magnetisation saturation for Fe and Fe/Pt layers are 
1700 x 10 3 and 690 x 10 3 A/m [19], respectively. 
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FIG. 1: Schematic representation of the position of the spins 
in the sample. 



In order to perform the dynamical calculations we need 
to obtain the equilibrium spin configuration, which is ob- 
tained by an energy minimisation routine without any 
static applied field. This is then taken as the initial con- 
dition for any subsequent dynamical calculations. The 
dynamics of interest in the Fe layer then arise due to the 
perturbation or excitation of the system upon the abrupt 
application of a magnetic field along some specific direc- 
tion with respect to the film plane. The evolution of the 
magnetisation is then evaluated as a function of time, 
from which we perform a Fast Fourier Transform (FFT) 
to obtain the relevant frequency spectrum. We have used 
the z-component of the magnetisation for all calculations 
of the spin dynamics since it is transversal to the direc- 
tion of the applied field, which is maintained in the x-y 
plane. 



III. RESULTS AND DISCUSSION 
A. Static configuration 

The simulations were performed using the conjugate 
gradient method in the OOMMF package to evaluate the 
spin orientation as a function of its position in the Fe 
layer. We have varied the Fe layer thickness (number of 
spins) in order to find an optimum thickness for which to 
obtain an exchange spring (ES) system (i.e. a rotational 
variation across the film thickness). The results for the 
rigid interface (first Fe spin fixed in the perpendicular 
direction) are shown in figure |2(a)| . Here we note that 



for thicknesses corresponding to 7 Fe spins, the system 
forms a rigid magnet (RM), whereby the Fe spins align 
along the FePt anisotropy (easy) axis as defined along 
the film normal. Adding one more spin we see that there 
is a significant relaxation to the spins in the Fe layer 
towards the in-plane direction. This relaxation further 
increases as further Fe spins are added. Almost full in- 
plane rotation (of the top Fe spin) occurs for around 25 
Fe spins. The transition from RM to ES is expected from 
energy considerations since at a certain thickness there 
are sufficient spins to allow the competition between the 
exchange coupling energy between the layers (causing an 
effective perpendicular anisotropy in the Fe layer) and 
the magnetostatic energy to be spread among the spins 
creating a domain wall. Full in-plane rotation can be seen 
as the construction of a 90° domain wall (DW). Such a 
transition has been predicted, for example by Asti et al. 
[6] and observed experimentally by Casoli et al. [12]. In 
addition to considering the rigid interface conditions, we 
have performed the corresponding calculations for a re- 
laxed interface. In this case we define the FePt layer to 
consist of 5 spins. This means that the domain wall can 
now penetrate the hard magnetic layer, helping to fur- 
ther reduce the tot al ene rgy of the system; the results 
are shown in figure |2(b)| Once again a transition from 
RM to ES is observed, though the number of Fe spins is 
reduced compared to the rigid interface condition. This 
can be expected since the extension of the DW into the 
FePt layer helps shift some of the energy from the Fe layer 
to that of the FePt. The static configuration is a sensi- 
tive function of the number of spins in both the Fe and 
FePt layers. We have performed extensive studies of such 
conditions and will be reported elsewhere [20 . For our 
purposes, we have chosen to study the Fe layer consisting 
of 25 spins using rigid interface conditions. This offers a 
system in which an almost complete 90° DW wall exists. 
We now proceed to study the magnetisation dynamics in 
this system. 



B. Magnetisation dynamics 

To study the magnetisation response to an abruptly 
applied magnetic field along some specific direction, we 
use the LLG, which is solved using the Runge-Kutta 
method, also in the OOMMF software. To understand 
the importance of material parameters and response of 
the system, we have made a series of simulations with 
applied magnetic fields of different strengths applied in 
various directions. The perpendicular anisotropy of the 
FePt layer provides the system with a unidirectional char- 
acter which will be transmitted to the Fe layer via the 
exchange coupling between the layers. Such properties 
will be expected to influence the dynamic properties of 
the Fe layer. This can be illustrated, for example, in the 
spectra obtained for the +20 and —20° directions, which 
are shown in figure [3(a)1 and[3(b7 



We note that the principle position (frequency) of the 
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(a) Angle with respect to the y direction as a function of the number of 
spins for a rigid interface. 
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(a)FFT for magnetic field of 1 T and 20°. 
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FIG. 2: (a) (b)|Show the spin angle as function of position in 



Fe layer. The first spin on each plot corresponds to the first 
unfixed spin in the chain. 



peaks in the spectra are the same; however, the relative 
intensities are very different - the principal peak in figure 
3(b) shows a higher intensity than in figure 3(a). There- 
fore this camouflages the importance of the other relevant 
peaks (frequency modes). Some of the principle charac- 
teristics of the dynamical properties can be elucidated 
from the study of the frequency field behaviour, which 
we have done as a funct ion of the direction of the applied 
magnetic field. In 4(a) and |4(b)j we show the dispersion 
relations of the dominant resonance frequency for various 
angles of the applied field. (We note that the 0° direction 
corresponds to the field applied in the plane of the film.) 
While the positive and negative directions appear to give 
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(b)Angle with respect to the y direction as a function of the number 
of spins for a 5 spins damped interface. 



FIG. 3: (a) (b) show the FFT spectrum of the Fe layer. 



very similar results, there are some important differences 
which are related to the unidirectional character of the 
system. The variation in each case appears as for a sys- 
tem with almost vanishing magnetic anisotropy, which is 
clearly not the case. 

In a system with a well defined anisotropy axis, the 
frequency field behaviour has a characteristic shape and 
depends on the type of magnetic anisotropy; uniaxial, cu- 
bic, etc. For example, in the case of uniaxial anisotropy, 
with the field applied along the hard axis, the disper- 
sion relation shows a marked cusp with a minimum cor- 
responding to the anisotropy field, Hk, of the system. 
While with the field applied along the easy axis the in- 
tercept with the frequency axis occurs to a resonance 
condition in which a field corresponding to the anisotropy 
field is applied; i.e. uj = jHk, see for example Vonsovskii 
[2T] , In the present case we see no well defined axis of 
anisotropy. This is because in our case this does not 
exist, since the directions of the spins in the Fe layer 
vary as a function of position across the thickness of the 
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(a) Plot of the resonance frequency as a function of the applied 
magnetic field for positive angles (measured in the first quadrant). 
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FIG. 4: (a) (b) depict the spectrum of frequencies as a func- 



tion of the applied magnetic field for several angles. 



film. This will mean that we cannot, strictly speaking, 
treat the system in the usual manner, despite the usual 
asym ptotic behaviour still being apparent [22 j. In figure 
5(a) we illustrate the resonance frequency as a function 
of angle (of the applied magnetic field) for 0.1, 1 and 2 
T. This shows a reasonably symmetric variation about 
the x-axis. However, a closer inspection reveals that for 
negative angles the frequencies are higher. 

A more sensitive display of this asymmetry is shown 
in the variation of the asymptotic slope o f the f requency 
field data as a function of angle, see figure [5(b)] The dis- 
parity between positive and negative directions becomes 
more marked for larger angles. In fact it would appear 
that there is a discontinuity in the variation in the region 
of 5°. This asymmetry and discontinuity, as well as the 
asymmetry in the angular dependence of the frequency 
are related to the unidirectional behaviour transmitted 
to the Fe layer from the FePt layer via the exchange in- 
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(a)Resonance frequency as a function of the angle for 
fixed applied magnetic field. 
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(b)Plot of the resonance frequency as a function of the applied 
magnetic field for negative angles (measured in the fourth quadrant). 




(b) Resonance frequency per applied magnetic field as a function of the 

angle. 



FIG. 5: (a)| (b) show the symmetry properties of the system. 



teraction. Added to this the rigid boundary conditions 
will be expected to accentuate this behaviour. 

To gain a fuller understanding of the dynamical re- 
sponse of the system to the abruptly applied fields in the 
simulation, it is instructive to consider the initial and fi- 
nal configurations; i.e. that of the initial state without an 
applied field (which will be the same for all simulations) 
and the final equilibrium state to which the system will 
relax. In figure [6] we show the spin configurations for spe- 
cific cases: the T case corresponds to the initial spin 
configuration used for each simulation. We also show 
the final state configurations for 1 and 2 T applied along 
the film plane (x-axis) and along the normal direction to 
the film plane (y-axis). In all cases the lower (Fe) spin 
is always in the perpendicular direction and arises due 
to the fixed (or rigid) boundary condition that we have 
imposed. In the case of a 2 T field applied along the 
film normal, all spins align along the field giving a sat- 
urated state. It will be noted that 1 T is not sufficient 
to push the DW out of the film and a small rotation is 
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(b) Variation of the angle between initial and final configurations 
as a function of the position of the spin in the sample in the y 
direction for 1 T and 2 T. 



FIG. 7: |(a) (b) show the relative deviations from the initial 
angle. The fluctuations observed in figure 7(a) are due to 
the error in determining the angle. Only the line of tendency 
should be considered and not the particular value of each 
point in the curve. 



FIG. 6: Initial spins configuration (0 T); final spin configu- 
ration in y (90°) and x (0°) direction. The first spin, in the 
botton of each chain, is the fixed spin from the FePt layer. 



still evident. Applying the field in the plane of the film 
effectively squeezes the 90° DW, reducing its thickness. 

The transition from initial to final state will define the 
dynamics of the magnetisation process via the preces- 
sional motion of each spin. There will be two main factors 
which will govern this motion; the effective field; H e ff, 
for each spin, which will depend on its position within the 
layer, and the overall angle difference between the initial 
and final states. The fact that the effective field varies 
as a function of position can be seen by the equilibrium 
spin distribution across the Fe layer. The difference be- 
tween initial and final state will govern in large part the 
amplitude of the motion of each spin in its movement; 
figure [7] shows these variations for the parallel and per- 
pendicular fields illustrated in figure [6] It will be noted 
that for in-plane fields the variation is negative, while 
for out of plane fields this will be positive (due to the 
convention we have chosen). This will therefore deter- 



mine the character of each spectrum and can be related 
to the resonance frequency and the relative intensities of 
the various peaks in the spectra. In addition to the am- 
plitude or precessional angle of motion, the intensity of a 
particular peak will also depend on the number of spins 
which are involved in that motion. To illustrate this in- 
terpretation we can consider some specific cases. Let us 
first consider the field applied in the plane of the sample 
(along the x-axis). The initial and final spin configura- 
tions are illustrated in figure [6| The first thing to note 
is that there is only a small change in the configuration 
itself, with the largest changes occurring at around the 
10th spin from the interface between Fe and FePt. For Fe 
layers of different thicknesses this will be different. We 
also note that the differences should increase with the 
strength of the applied field; see figure [7(a)] When the 
magnetic field is applied along the film normal we notice 
that the character of the spectrum becomes much more 
regular; the spectrum for a 2 T field applied in the per- 
pendicular direction is shown in figure [8] Similar spectra 
are obtained for other field values. In this case the dif- 
ferences from initial to final state are much greater, and 
increase as we move from the interfa ce to the outer sur- 
face of the Fe layer, as shown in figure (7(b)) In traditional 
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FIG. 8: Spectrum of frequencies obtained by FFT for an ap- 
plied magnetic field of 2 T and 90 (with respect to the y 
direction). 



magnetodynamic measurements performed by ferromag- 
netic resonance (FMR), for example, multi-peaked spec- 
tra are usually associated with the existence of standing 
spin wave resonance modes which arise from the mag- 
netic confinement in the perpendicular direction and the 
specific boundary conditions which apply. Such consid- 
erations lead to a resonance equation for the resonance 
frequency in function the applied magnetic field and the 
spin wave wave- vector of the form of the Kittel equation 
[23] - see equation [5] 

If we consider that the spins at the lower interface are 
fixed in the perpendicular direction, while those on the 
outer surface are free to move (as in the bulk), we can 
easily arrive at the allowed values for the spin wave modes 
as given by: 



(2n^l) 
k ~ 2L W 



n = 1, 2, 3... , 



(7) 



where n denotes the mode number and L the thickness 
of the Fe layer. In this case both odd and even modes 
should be allowed since the boundary conditions are non- 
symmetric. A plot of frequency against the square of the 
wave- vector should yield a linear variation - according to 
equation [5] Using our data we have made such plots for 
various applied fields along the perpendicular direction; 
these are shown in figure [9] We note that while much of 
the data fit on linear portions, the extremal values are 
clearly not in agreement. Therefore we believe that the 
various peaks in the spectra do not arise from the excita- 
tion of standing spin wave modes. Rather, we think that 
the existence of the multi-peaked spectra is due to the 
different local conditions (local effective fields) of each 
spin, there will be a distribution of resonance conditions 
which are sufficiently spread out as to be distinguishable 
resonances. These local effective fields will reflect the 
differing conditions and the local environment, where a 



local energy minimum is encountered, giving rise to the 
direction of the spin at that position. The intensity dis- 
tribution also reflects the initial to final state directions, 
which in much of the layer are quite large (see figures 
|6]and|5|. This will mean that the initial precessional 
angles will be large and hence give increased intensities, 
which are proportional to the transversal component of 
the magnetisation. 
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IV. CONCLUSIONS 

We have made a detailed study of the model system 
comprised of a hard magnetic layer (FePt) with perpen- 
dicular anisotropy exchange coupled to a soft magnetic 
layer (Fe). The first step to the study of the magneti- 
sation dynamics in this system concerns the determina- 
tion of the equilibrium condition. This is achieved by 
minimising the free energy density of the system and ob- 
taining the orientation of the individual local spin ori- 
entations. We have performed such calculations for the 
system with varying Fe thicknesses. Since we are inter- 
ested in studying the exchange spring system, we use 
these calculations to choose the Fe thickness we require; 
i.e. where there is significant rotation of the magnetic 
moment through the film. We see the transition from 
the RM to ES occurs for a system with a thickness cor- 
responding to 6 and 7 atomic spins for rigid and free 
boundary conditions, respectively. This corresponds to a 
thickness of around 0.86 and 1.0 nm and compares rea- 
sonably well with the calculations of Asti et al. [6 , who 
obtain a value of between 0.6 and 0.8 nm, depending on 
the thickness of the FePt layer. 

The dynamics of the system - in response to an 
abruptly applied magnetic field - were studied as a func- 
tion of the strength and direction of the applied field. 
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This allowed us to obtain the frequency spectrum and 
therefore construct the dispersion relation. We have con- 
sidered only the dynamics of the Fe layer, since this 
should be separate from the FePt dynamic response due 
to its elevated anisotropy, and is supported by FMR mea- 
surements which only detect Fe resonance features at low 
frequencies [13]. We note that the lack of a well defined 
anisotropy axis in the Fe layer means that the angular 
variation of the frequency field characteristic differ from 
that of a normal film, and always pass through the origin. 
The high field asymptotic gradient varies with direction, 
whose angular dependence is related to the exchange cou- 
pling with the hard layer. A small asymmetry is noted 
around the x-axis which is due to the unidirectional char- 
acter of this exchange coupling. For magnetic field ap- 



plied in the perpendicular direction we obtain a rich fre- 
quency spectrum which can be understood in terms of 
the variation of the local effective magnetic field across 
the Fe layer. 

As an aside, it is worth noting that the OOMMF soft- 
ware provides sufficient resolution to allow us to treat the 
system as atomic layers of spins in the direction perpen- 
dicular to the film plane. This is supported, for example, 
by the fact that we were able to reproduce the RM-ES 
transition, and is supported by experimental measure- 
ments [12] and other calculations [20] . 

Acknowledgements. We thank OOMMF mailing 
list for insightful discussions. 
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